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Two light pulses propagating with ultra-slow group velocities in a coher- 
ently prepared atomic gas exhibit dissipation-free nonlinear coupling of an 
unprecedented strength. This enables a single-photon pulse to coherently 
control or manipulate the quantum state of the other. Processes of this kind 



\^ . result in generation of entangled states of radiation field and open up new 

prospectives for quantum information processing. 
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It has been known for more than thirty years that hght fields or photons can interact with 
each other in atomic media much hke massive particles do |1|]. However, the strength of the 
interaction of two single light quanta is typically extremely weak. As a result, conventional 
nonlinear optics is feasible only when powerful laser beams, containing a large number of 
photons, interact in nonlinear materials. 

This Letter describes a method that allows for two slow light pulses |@-§] of tiny en- 
ergies to interact in a resonant ensemble of atoms. When Electromagnetically Induced 
Transparency (EIT) [^^j is established for both of these pulses, they will propagate with 
ultra-slow but equal group velocities, and a very efficient nonlinear interaction between 
them will take place [0. This interaction can be maintained for a very long time without 
dissipation, resulting in a new regime of quantum nonlinear optics. Specifically we describe 
here a scheme in which a traveling light pulse with an energy corresponding to that of a 
"single photon" {hu) can modify the refractive index of a second traveling pulse such that 
the latter experiences a non-linear phase shift on the order of vr. This process allows one to 
create strongly correlated (i.e. entangled) states ^ of interacting photons and to generate 
macroscopic quantum superpositions such as Schrodinger cat states |g]. The technique de- 
scribed here opens up interesting prospects for coherent processing of quantum information 
lOHTTl in applications such as quantum computation and quantum communication. 



It has been widely accepted that the use of ultra-high finesse micro-cavities is essential 



for having strong interactions between single photons [jT2|-[Tj]. In the present contribution 
a radically different approach is taken, in that the properties of an optical material are 
"designed" to achieve an unprecedented strength of nonlinearity without the need for a 
cavity. Earlier work has already demonstrated that such a coherently controlled dense atomic 
ensemble displays remarkable phenomena [§],§[, which are of great interest for nonlinear 
optics at a low light level |T5|-p!^. 

We consider a process in which a pair of very weak optical fields Ei, E2 interact resonantly 
with an ensemble of A^i multi-state atoms (species A) as depicted in Fig.l. These atoms 
are coherently driven by a classical (i.e. many-photon) light field with a Rabi frequency Qi 



tuned to resonance with an atomic transition |ci) -^ |ai). We assume that the frequency 
of one of the weak optical fields {Ei) is tuned to resonance with transition \bi) -^ |ai). 
Quantum interference induced by the driving field results in a sharp transmission resonance 
in the spectrum of the weak field Ei (Fig. 2). The second weak field {E2) couples another 
optically allowed transition |ci) -^ \di) with a single photon detuning A. When E2 is 
absent a properly tuned Ei propagates without loss or refraction, but its group velocity is 
substantially reduced: this is the essence of EIT. The off-resonant field E2 induces a Stark 
shift of the state |ci) and therefore modifies the refractive index of the first field. Resulting 
phenomenon corresponds to cross-phase modulation of two weak light waves : since the 
dispersion of refractive index is very steep (Fig. 2), small Stark shifts result in a large index 
change. 

The key idea of the present work is to arrange the conditions such that the two weak 
light pulses can interact in transparent, nonlinear media for a very long time. This is 
the case when pulses propagate with equal, slow group velocities. Photon interaction will 
then cause large nonlinear phase shifts even if pulses of ultra-small energies are involved. 
A possible technique to achieve equal group velocities for two weak pulses is explained in 
Fig.l: here N2 A-type atoms of a second kind (B) are added to the interaction region, 
and the frequency of the field E2 is chosen to be on resonance with the optical transition 
1^2) -^ \0'2)- Application of a second driving classical field ^2 of appropriate intensity will 
result in induced transparency and properly reduced group velocity for the field E2. The 
mixture of atoms required to achieve such group velocity matching can be designed by using 
different isotopes of alkali atoms, and by applying an appropriate magnetic field to create 
the desired level splittings. 

We describe quantum radiation by multi-mode field operators 0: ^12(2:, t) = 



2Trc/LJ2k(^ki,2{'t)e~^^^, where aki is the annihilation operator corresponding to the field 
mode with wave- vector k. L is the quantization length and c is the speed of light in vacuum. 
In the present approach we choose to restrict the quantization to a continuum of modes with 
wave-vectors ranging from k^ — Auj/2c to k^ + Auj/2c, which in turn gives equal space-time 



commutation relations [Ei{z,t),Ej~{z,t)] = Auj6ij, that depend explicitly on the bandwidth 
of the system Auj. This bandwidth can be determined, e.g., by the detection process. As 
we will discuss shortly, the spectral width of the transparency window also gives a natural 
limitation on the maximal bandwidth. 

The dynamics of the atomic media is described by Heisenberg equations for atomic 
polarizations and coherences. The finite bandwidth of quantized fields allows us to apply 
adiabatic elimination of the atomic degrees of freedom. By disregarding time derivatives of 
third and higher order, we arrive at the following evolution equations for the field operators: 

+ iriE+E2Ei + Fi (1) 

+ ir]EtEiE2 + h. (2) 



The adiabatic expansion procedure used in derivation of Eqs. (Q,^ is analogous to that 
of Refs. pr5| , pT| for classical fields. Here, Vg = c/(l + Ug) is the group velocity for the 
corresponding pulse, Ug = g'^N/\Q\'^, and n being the rate of residual single-photon loss 
K = Ug'jhc/c. Fi and F2 are delta-correlated noise operators associated with dissipation. 
T] = [rigg'^ / {i'jcd + A)]//(27rc2) is the rate of nonhnear interaction between pulses and P = 



^g7afe/(|^Pc). g = J^abO'c/{2Al) IS a normalized atom-field coupling constant, A is the 
cross-sectional area of the quantized fields, a is resonant absorption cross-section and / is 
the length of the interaction region, 'jij are linewidths of \i) -^ \j) atomic transitions. We 
assumed that all interacting atoms have identical linewidths and coupling constants, and 
that the number of atoms are equal {Ni = N2 = N). To obtain equal group velocities, 
identical Rabi- frequencies for the driving fields were chosen {Qi = Q2 = ^) ■ 

Eqs. (|1],0) together with the commutation relation define a quantum field theory of 
interacting pulses in a coherently prepared medium. In cases when losses can be neglected 
this field theory is analogous to that of a weakly interacting Bose-gas, but with an imaginary 
effective mass term. This imaginary mass term results from the finite bandwidth of the 



transparency window and causes a spreading of the input pulses (Fig. 2). 

We first consider the classical limit of this theory by replacing operators Ei^2 with their 
expectation values -Ei,2- In the ideal case A can be chosen to be large and hence rj is purely 
real. When attenuation and pulse spreading are small enough to be neglected, Eqs. (^0) 
can be easily solved by Ei^2{z,t) = i?i_2(0, t') exp(2?7|_E'2,i(0, t')pz), which corresponds to 
cross-phase modulation of the weak beams. Note that this result is expressed in terms of 
retarded time t' = t — z/vg. Magnitude of the nonlinear phase shift at the peak of the pulse 
can be expressed in terms of energy E and duration T of a Gaussian pulse: 



phase shift = ./^^^#I^. (3) 

Here nonlinearity rj is expressed in terms of group delay Tg = Ugz/c. In the limit when 
Vg ^ c the later corresponds to the interaction time of the pulse with the medium. 

It is apparent that large phase shifts exceeding vr at energies corresponding to a fraction 
of hv appear to be possible when the interaction time exceeds the pulse duration of a 
tightly focused (cr ~ A) laser beam. Before proceeding it is important to re-examine the 
assumptions resulting in this striking conclusion. 

First, we note that a stringent limitation imposed by the difference of group velocities of 
interacting pulses Jl^] is absent in our case since the interaction is designed such that group 



velocities are equal. Second, even in the absence of various absorption mechanisms, nonlinear 
phase shifts are limited by the bandwidth of the transparency window which decreases with 
propagation distance |l^ AoJmax = {13 z)'^/"^. After a sufficiently long propagation this 
results in spreading of the pulse (see Fig.2). In order to avoid losses due to spreading, Tg/T 
should be smaller than 



9?T„ aN ,^. 



The quantity under square root corresponds to an optical depth of the medium |]19| . Hence 
it is essential to operate in optically dense gas. Finally we note that the interaction time (r^) 
is itself limited due to absorption associated with decoherence of the "dark state" 0; i.e. 



Tg < 7^^ must be satisfied. This implies that EIT resonances should be strongly saturated, 
i.e. ^"^ I'^ahlbc > Na /A ^ 1, in order for nonlinear phase shifts to be large. 

We now consider quantum dynamics of the nonlinear interaction. When absorption is 
negligible and the bandwidth 5 of the pulses satisfy 5 < Au; < ^uj^ax^ quantum Eqs. (^0) 
can be solved by 

^1,2(2;, t) = E,^2{t') exp [i 7] El,{t')E2,i{t')l] , (5) 

where Ei{t) are Heisenberg operators describing the input fields at z = 0. We first analyze 
the evolution of most "classical" of all possible input wave packets. These are multi-mode 
coherent states P \ai, 02) = Y[k,j |q^i) ^ |q;2)) which are eigenstates of input operators Ei{t) 



with eigenvalues (at z = 0) ai{t) = J2Trc/LJ2k'^i^^''^^- After propagation through the 
nonlinear dispersive medium, the following expectation values of the fields can be measured: 



{Ei,2{z,t)) =ai,2(t') X 



exp 



-2sin2($/2)+isin($)] l"2,i(OI 



(6) 



Aa; 

where quantum phase shift $ = [l/(47r)cr/y47c(i/A]Aa;rg. These solutions have the same 
form as those obtained in Ref. |^ for single-mode fields and generalize that earlier result 
to the multi-mode case that is appropriate for the traveling-wave geometry. 

Eq. (P) reproduces the classical result only when $ — * 0; for sufficiently large interaction 
times the quantum dynamics of wave packets deviates substantially from the classical case. 
In particular, both phases and amplitudes given by Eq. (|]) exhibit periodic collapses and 
revivals. The origin of this behavior can be understood by noting that each component of 
input coherent states |ai)|a2) is itself a coherent superposition of many Fock components 
In^rrag). During nonlinear interaction each of these components acquires a different phase 
change resulting, for sufficiently large $, in quantum dephasing of the original coherent 
states. However when $ reaches multiples of 2tt, all components "re-phase" such that 
original coherent state is reproduced. At intermediate values of interaction times, re-phasing 
to other macroscopic states can occur. For example, when $ = vr the output state of two 
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fields can be verified to be: I?/') = | (|ai, 02) + | — ai, 02) + |«i, — ^2) — I — cti; — «2))- This 
is an entangled superposition of macroscopically distinguishable states. Superpositions of 
this kind have no classical counter-part, and correspond to Schrodinger cat-like states [j^. 

Consider now a different type of input quantum state corresponding to multi-mode single- 
photon wave-packets: |lj) = J2k^kO'ki\0) where Fourier amplitudes ^k are normalized such 
that J2k l^fcP = 1- In free space these single-photon states represent traveling waves initially 
(t = 0) localized around z = 0. The dynamics of individual wavepackets is fully described by 
a single photon "wave-function" '^i{t,z) = {0\Ei{t)\li) . Correlations of these wavepackets 
emerging due to photon-photon interaction are described, in turn, by correlation amplitudes 
^uiti, Zi;t2, Z2) = (0|^i(ti, ^1)^2(^2, -22)|lil2)- After propagation though the atomic cell 
the correlation amplitude is given by: 

^u{t';,z, ■ t'2',^2) = ^i(o,t';)^2(o,4') 

X "1 + sinc[Aio{t'; - t2)/2](e^* - 1)] , (7) 

where t'- = ti — l/vg — {zi — l)/c. Eq. (|^) indicates that the nonlinear interaction alters the 
mode structure of the pulse in addition to generating a phase shift on existing components. 
This phenomenon is related to the classical effect of pulse broadening due to cross-phase 
modulation. It is important however that equal-time ("coincidence") correlations (t'/ — 
t2)A(jj <^ 1 indicate that a pair of single photons acquire a phase shift $ as a result of the 
nonlinear interaction. The magnitude of $ can easily exceed tt, when TgAumax ^ 1- 

These large nonlinear phase shifts can be used to create quantum entanglement. For ex- 
ample, if the photons are initially in coherent superpositions of two states (for example the 
two polarization states |1±)), and only one of these states (|1+)) is subject to the strong non- 
linear interaction, the resulting state cannot be factorized into a product state of individual 



modes [£^ . Operations of this kind form the essence of quantum information processing |]T0 
The example considered above indicates, however, that the present approach to quantum 
entanglement differs conceptually from the techniques discussed previously. In particular, 
usual approaches to quantum processing are based on systems whose Hilbert space can be 



restricted to a 2 dimensional subspace (qubit). In contrast, the present technique involves 
multi-mode, traveling wave excitations where the 2-dimensional Hilbert space of the qubit 
(i.e. polarization) and the external degrees of freedom of the photon field become coupled. 
It is interesting to consider how various concepts of quantum information theory can be 
applied in the context of our approach. 

It is clear that the feasibility of quantum entanglement depends upon the large delay- 
bandwidth products l\ijJraaxTg Corresponding to ultra-slow pulses. In experiments involving 
ultra-cold |@| and hot |l^ atoms values of this product on the order of few tens have been 



observed. Potentially, an increase by two to three orders of magnitude is likely, which should 
allow for a high fidelity of entanglement. We conclude by noting that exceptionally large 
nonlinearities have already been measured in various experiments P,|3|,p!9| . In particular, an 
efficient nonlinear phenomenon corresponding to vr phase shift has been observed with pulse 



energies corresponding to less than 10 photons per atomic cross section ||22[ . 
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FIGURES 

FIG. 1. Prototype composite atomic medium for strong nonlinear interactions at the sin- 
gle-photon level. By adjusting detunings such that the two weak fields Ei,E2 and two classical 
fields r2i,r22, are in resonance with transitions in atoms A and B as shown, resonant absorption 
can be eliminated and the group velocities of two weak pulses can be made equal. The scheme 
can be implemented, for example, using a natural mixture of atomic rubidium. It consists of two 
isotopes Rb^'^ (natural abundance ~ 25%) and Rb^^ {^ 75%) with transition frequencies differing 
due to isotope shifts and nuclear spins. Propagation of tightly focused beams over long distance can 
be achieved using non-diffracting beams of Bessel shape, fiber waveguides, or focusing properties 
associated with FIT to induce waveguiding optically. 

FIG. 2. (a) Spectrum of transmission and refractive index corresponding to FIT. Quantum 
interference in the atomic ensemble induced by a coherent driving field creates a sharp resonance 
in the transmission of a weak optical field, accompanied by rapid variation of the refractive index 
(red curve). This rapid variation causes a dramatic reduction of group velocity. The presence of 
a second weak field causes an effective shift of the resonant frequency (green curve) , which results 
in corresponding change of a refractive index and hence the phase of the first weak field, (b) 
Propagation dynamics in coherent media: when resonant pulses enter the medium, they exhibit 
spatial compression and very soon slip behind the reference pulse (white) that does not interact 
with the atoms. The total time that two slow pulses can spend in the medium is limited by the 
residual single photon loss, and by the spreading of the pulses after r^ ~ i^'^T'^/jab- The spreading 
is due to a finite bandwidth Aujrnax of the FIT resonance. 
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